Quantum walk on a line for a trapped ion 
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We show that a multi-step quantum walk can be realized for a single trapped ion with interpolation 
between quantum and random walk achieved by randomizing the generalized Hadamard coin flip 
phase. The signature of the quantum walk is manifested not only in the ion's position but also its 
phonon number, which makes an ion trap implementation of the quantum walk feasible. 

PACS numbers: 42.50.Ex,03.67.Ac,03.65.Yz 



Introduction:- Quantum information proeessing 
promises revolutionary advances in communication 
and computing with secure long-distance quantum key 
distribution [l| and quantum computing [4] as two im- 
portant long-term goals. In the medium term, progress 
in quantum information implementations has been most 
pronounced for quantum communication protocols 
which consume entanglement to enable quantum- 
enhanced communication. An experimental quantum 
walk (QW) 0, i, i would be a major advance forward 
towards the ultimate goals of quantum information 
processing: benchmarking coupling between qubit and 
'bus' mode, test models of the environment and control- 
ling decoherence, and simulate exponentially-enhanced 
quantum algorithms over classical counterparts Our 
aim is to realize a multi-step coined QW, which would 
be implemented in an ion trap Q . 

The random walk (RW) is ubiquitous in physics, chem- 
istry, mathematics, and computer science: this process 
underpins Brownian motion and diffusion process, is used 
in satisfiability proofs (SAT), and is intimately connected 
with the Wiener measure Q . Quantization of the RW Q 
has led to new quantum algorithms Q and fascinat- 
ing physics such as decoherence-induced diffusion reduc- 
tion [6[ . Our goal is to see the QW realized in the labo- 
ratory. For this purpose, we consider the simplest RW: a 
single walker whose two degrees of freedom are position 
on regular one-dimensional lattice and a single two-sided 
coin that generates random bits. Each coin flip generates 
a result or 1, causing the walker to step left or right 
on the line, respectively. In the quantum version of this 
walk on a line, each position is a state in Hilbert space, 
and the coin is a qubit whose flip is a unitary evolution: 
superpositions of position and entanglement between the 
coin state and the walker's position are now possible. In 
the QW, the walker commences at a point on the line 
and alternates between unitary coin flipping and making 
left or right steps that are entangled with the coin state. 

This single-walker QW has not yet been realized ex- 
perimentally. Most proposals focus on implementing the 
QW in phase space (walking around a circle in an ab- 
stract position- momentum space) Q , and even this easier 
case (due to wandering in position space being strongly 



bounded) has not been realized. Here we show that the 
QW on a line is indeed achievable in an ion trap. 

Travaglione and Milburn (TM) first proposed a QW 
implementation: the walker's degree of freedom would 
be the position state |a;) of a single trapped ion and 
the coin state corresponds to the up state |t) and down 
state II) of the ion's electronic degree of freedom. Al- 
though a seminal proposal, unfortunately it is not vi- 
able, due mainly to four drawbacks: (i) impracticality 
of measuring x; (it) unavoidability of higher-order Lamb- 
Dicke (LD) contributions for large numbers of steps [l^l ; 
(iii) the need to control decoherence to enable interpo- 
lation between the RW and the QW so that comple- 
mentarity, hence quantumness, can be verified and 
most importantly (iv) the impossibility of reading the 
coin state (projection onto \ i) or |t)), because this proce- 
dure requires scattering many photons on the ion, which 
inevitably destroys its motional state. One suggestion for 
addressing the impracticality of a;-measurement by using 
instead a quantum network and multiple ions is inter- 
esting but also impractical, and suggestions of quantum 
walks on circles in phase space are also of value but avoid 
entirely the core issue of realizing a QW on a line [5|, i6|] . 

TM's concept of an ion trap implementation of the 
QW is laudable, but clearly major advances are required 
to bring their concept of experimentally realizing a QW 
on a line to feasibility. Here we overcome all four draw- 
backs: (i) we replace x-measurement by measuring in- 
stead phonon number n via Rabi flopping and show that 
the phonon counts exhibit an unambiguous QW signa- 
ture; (ii) we include higher-order LD contributions and 
show how they lead to breakdown of the QW on the 
line; (iii) we introduce a random-phase generalization to 
the coin flip and show how averaging enables interpola- 
tion between the QW and the RW [ll|, and (iv) devise 
a method for inferring the phonon number distribution 
from electron shelving by looking at the carrier transi- 
tion instead of the first blue sideband [l^l- In making 
the advance from concept to design of an ion-trap imple- 
mentation for a QW, we introduce some new and valuable 
methods that have broad beneflcial implications for ion- 
trap-based quantum information processing in general. 

Ideal QW on a line:- Before discussing the full, feasi- 
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FIG. 1: (a) Mean phonon number and variances of posi- 
tion and momentum with a — 0.565 up to A'^ = 17 steps, 
(b) Wigner function W{x,p) for tire walker state at = 10. 



ble implementation of the QW on a line in an ion trap, 
here we study the ideal QW on a line and show that 
the phonon number distribution P„ carries the signature 
of the QW as well as the inaccessible position distribu- 
tion P{x). 

Sufficient criteria for experimentally demonstrating the 
QW arc: a single walker whose position x is restricted to 
a one-dimensional lattice € Z with step size a between 
lattice points; the position is incremented k k + l if the 
coin state is |t) and decremented k i-^ k—1 for capa- 
bility of unitarily flipping the coin such that ||), | J,) evolve 
to equal superpositions thereof; quadratic enhancement 
of spreading for the QW compared to the RW; and con- 
trollable decohcrence to interpolate between the QW and 
the RW. 

Mathematically the unitary coin toss operator C is 
given by •\/2C'(0) = l — iax cos (p+iay sin for 1 the iden- 
tity and ffx (flip), (phase-flip), and ct^ = diag(l,— 1) 
(phase gate) designating Pauli operators. The phase (j) is 
arbitrary but must be constant for a unitarily evolving 
QW. The walker's step to the left or right, entangled by 
the coin state, is enforced by the unitary operator T = 
exp (iap (8) (Tz) for p the position translation generator 
(i.e., a momentum operator) and a the step size in a line. 
Each step is effected by (5(0) = T [1 (g) C(<^)], and the 
evolution to N steps is {{M) = U'Li T\t(» C((/)f)]. 
For the ideal QW, (fig is constant over all steps and typi- 
cally assigned a value of 7r/2. For ak = {N — 2k)a, an ini- 
tial walker-|-coin state l-^o) = |0) (ID -(- It))/ -v/S evolves to 

li'N) = Q^IV-o) 
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walker state is pN = Tr (|?/'jv) (V'Jvl), which has position 
distribution P{x) = {x\ pn \x) and phonon number dis- 
tribution Pn = (n\ \n). Position variance cr^(x), mo- 
mentum variance cr^ (p) , and mean phonon number n are 
shown in Fig. [ija): evidently cr^(x) oc N"^ and cr^(p) ~ 
constant up to iV = 17 steps; for the RW, cr^(x) oc N, 
and this quadratic enhancement of position spreading is 
a signature of the QW. 

Fig.[T]Ja) also shows that the mean phonon number h on 
N"^ , which we identify as an alternative QW signature: 



our choice of initial state yields n = [cr^(x) -I- cr^(p)] /2 a 
iV^, whereas, for the RW, n (x N . Therefore, a quadratic 
enhancement in spreading of n is just as good a signature 
of the QW as is enhanced spreading of a'^{x)\ Fig. [Hb) 
reveals the invariance of the state under x —x or p i— > 
—p mappings by depicting the Wigner function |14{ for 
the state after ten steps. Position x and momentum p 
distributions are marginal distributions of W{x^p) 

Controlling decohcrence:- Although decohcrence oc- 
curs naturally in the laboratory, for example due to 
magnetic field fluctuations, which generate random z- 
rotations, this undesirable decoherence is largely elimi- 
nated by spin echos (built into our technique). To intro- 
duce controllable decoherence, we uniformly randomly 
choose each G (— Tr/g, Tr/q) with q > 1 a. controllable 
parameter that yields the RW for q = 1 and the QW 
for q oo. The resultant walker state is p{{(j)e}). De- 
coherence is achieved by running the experiment many 
times with different random sequences i{(j)e}) each time, 
then averaging the density matrix at the N*"^ step of 
each run to obtain ppj. For a(x) the position-spread 
and n the mean phonon number for pTv, we conjecture 
the power-law rules lna(a;) cx <;lnA^ and Inn cx ^InA'^, 
with (<; w 1/2, C w 1) for g = 1 and (<^ ^ 1,^ ^ 2) for 
q oo. Thus, q controls decoherence and interpolates 
between the two extremes of QW and RW. 

Ion trap implementation:- A single trapped ion (e.g., 
^Be"'') is confined in a radio frequency (RF) ion trap. 
Electronic (coin) and motional (walker) degrees of free- 
dom are coupled by 'carrier' and 'displacement' laser 
beams [15[. The carrier-beam difference frequency is set 
to the frequency difference of the coin states {|j) , |t)}. 
The difference frequency S of the 'displacement' Raman 
beams is close to the ion's motional mode frequency uj^. 

The initial state is prepared by laser-cooling the ion 
to the motional and electronic ground state, |0) ||), 
then applying a 7r/2 pulse, which creates an equal 
superposition of \l) and |t). By applying displace- 
ment Raman beams, the interaction Hamiltonian is 

H, = [c-H**-^)i^(i,ye-^*) + he] (n^ li) (il + IT) (Tl 

for D{a) = exp(aa^ — a*a) the unitary displacement op- 
erator and the carrier Rabi frequencies f2| = — ri|/2: 
This interaction approximates the desired evolution. 

We expand Hi perturbatively in powers of the LD pa- 
rameter r]] this expansion is valid provided that ?7 |(a) | 
and I (a^) | — |(a)|^ = (j'^{x) + (t'^{p) + 2Gov{x,p) are small, 
so (T^(a;), (7^(p) and Cov(a;,p) are each small. Thus, 
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Evolution U = cxp 
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over time t is ap- 
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proximately 

D{2n^r]t)Bi2n^r]h)UoB{2n^) ID (il 

+D{-n^r^t)B{-n^7jh)Uoff{~n^) It) (Tl + oiv^), (2) 

with each product of unitary operators in the sum com- 
prising resonant unitary evolutions, whose exponents are 
linear in t, and non-rcsonant evolution, with rapid terms 
such as exp (iu>zt). 

Evolution ([2]) is dominated by resonant terms 
of the Taylor expansion of Hi such as B{P) = 

g/3/6[(at)^a+ataat+o(at)2]_hc ^ g/3/2at -hc^ which 

arises from the 3'''^-order resonant terms, and 
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to 0{r]'^). Here S{z) = exp |iz*a^ — he}, for z ~ 
, is the squeezing operator, which 



arises from 2"'^-order resonant terms of the Taylor ex- 
pansion for Compared to resonant terms, the effect 
of C/off on n and cr(x) is small. We also neglect commuta- 
tors arising from expanding exp |— i i/i(t')dt'| except 

ioi D, B and S, which are non-negligible. We now clearly 
understand the small and large contributions to U. 

Without B and Uos, Eq. ^ is essentially a displace- 
ment of 2f2|77i or ^Q^rjt if the coin is | or "f , respectively. 
These asymmetric steps can be replaced by identical left- 
ward and rightward steps by alternating two 7r-pulses on 
the spins with two displacement steps U and W to yield 
a desirable evolution from U: Utot = (1 (8) X)W{1 (g) 
X)U{t ® C), with X = \i) (Tl + IT) (T|. Ignoring higher- 
order on-resonant term B and off-resonant term UoS 
yields Utot « D{3n^7jt) |i)(i| + D{-3n^r,t) |T)(T| to ob- 
tain a symmetric step size of ~ i^n^rjt. 

In Fig. ^a) we see that the position peaks re- 
main closer to the center compared to the ideal QW. 
The change in displacement of peaks is small but 
nonnegligible as we can see from the following argu- 
ment: The evolution is dominated by displacements 
D{3ft^r]t) K 1 + 3r2f?7i(a'l' — a), with first-order term 
~ 3r2|7?iv^ for I (a) I = = 3A^rif77t, whereas the S''*- 
order B{±3Vl1r]^^t) contribution scales as ilit{r]y/E)^ /2, 
which is responsible for fi-dependent displacement. Thus, 
the LD parameter needs to be kept small to ensure the 
largest number of possible steps iVmax with quadratic en- 



hancement of spreading: 77 ^ y ^/ y/ N^i^x^-^t. 

Whereas higher-order resonant terms modify the QW 
as shown in Fig. [2fa), non-resonant effects are much 
smaller than resonant contributions. We test the contri- 
bution of each unitary operator in UoS by forcing some 
unitary operators to be identities; then we evaluate how 
much each of those operators affects the dynamics. By 
this procedure we ascertain that non-resonant opera- 
tor contributions D ^^ifl-^r]^^^^ and B (^iifli^rj'^^^p^ 




FIG. 2: Quantum walker after A'^ = 10 steps, (a) Position 
distribution for the ideal QW on a line (black) with a = 
0.565('-^ 3Q,^r]t) and for the ion trap implementation (red) 
with parameters {5,u!z,^] )/2Tr — (4, 4, 0.3)MHz, 77 = 0.1 and 
pulse duration t = Ifia. (b) Contour plot of Wigner function 
(min=— 0.095, max=0.052) for ion-trap implementation, with 
each contour corresponding to a step of 0.025. 



arc primarily responsible for creating the pincer-like mo- 
mentum sidebands observed in the walker's Wigner func- 
tion shown in Fig.[2{b). These momentum sidebands are 
directly responsible for the small monotonic increase of 
momentum variance observed in Fig. [Ija), whereas mo- 
mentum variance is constant for the ideal QW on a line. 

We establish numerically that squeezing S in the evo- 
lution is responsible for slight asymmetry of the position 
distribution in Fig. [IJa) . This asymmetry arises because 
the position distribution peaks are squeezed conditioned 
on the coin state; when probability amplitudes are co- 
herently added, the position distribution symmetry is 
slightly violated. Off-resonant contributions also effect 
a small rotation of the Wigner function in phase space 
due to a e'^" contribution to Utot for small 9. 

Counting phonons:- Previously [Tsj the motional num- 
ber distribution has been determined by driving the ion 
on the first blue sideband and Fourier transforming the 
atomic population in the |T) as a function of drive dura- 



tion Px(t) = 1/2 1 + Y.n=i, Pn cos(f) 



n,n+l 



for Pn the 

> IT,« + i) 



n-phonon probability and ^n.n+i the |T,n) 
Rabi frequency. 

Outside the LD regime, ^n,n+i "y^f^nl'?^); where 
the generalized Laguerre polynomial is nonmono- 

tonic in n thereby leading to ambiguities in determin- 
ing Pn by this method. This problem can be redressed 
for low 7] by Fourier transforming the carrier signal for 
the transition \l,n) ^ \^,n)) where 17„^„ oc i°(7y^) [isj . 
These frequencies fln,n are monotonous and distinguish- 
able for n < 60, which is a sufficient range to observe the 
hallmarks of QWs vs RWs. For 7] < 0.2, the carrier alone 
is sufficient to find P„ for n < 25. However, once P„ is 
known for n < 25 this information can serve to lift the 
ambiguities on the blue sideband, which then in turn can 
be used to determine P„ for 25 < 71 < 60. 

Simulations:- We simulate the walker -|-coin dynamics 
and calculate the mean phonon number n and position 
spread ^{x) in two cases: the ideal walk on a line in 
Figs.[31Ja, b) and the walker as a trapped ion in Figs.[3Jc, 
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FIG. 3: (color online) In-ln plot of ideal walker's (a) mean 
phonon number n and (b) position spread a-{x) and ion trap 
realization of walker's (c) mean phonon number and (d) po- 
sition spread as a function of step number N for decoherence 
q ~ 1 (green), q = 5 (blue), and g = 20 (red). Experimental 
parameters and a for the ion trap realization (c, d) are as in 
Fig. 2. The slopes of the dotted lines, which interpolate from 
the ideal QW to RW (shown in each figure as solid lines with 
greatest and least slope, respectively), are (a) 1.003, 1.660, 
and 1.938, (b) 0.510, 0.860, 0.990, (c) 1.061, 1.530, and 1.902, 
and (d) 0.496, 0.824, and 0.985. 

d). Decoherence is controlled by q: RW for g — > 1; QW 
for q oo. Step size a = 0.565 corresponds to 3f2|ryt. 

The In-ln plots reveal the small- and large-q power law 
relationship between either n or a{x) and step number TV, 
as predicted; furthermore the slopes approach 1 for the 
RW and 2 for the QW, thus confirming that both n and 
^(x) suffice for observing the RW-QW transition through 
controlling q. Slopes of h and d'{x) for the ion trap case 
are slightly smaller than for the ideal QW, but this small 
degradation is fully explained, namely nonlinear and non- 
resonant contributions to the evolution. The RW-QW 
transition is excellent despite these pragmatic considera- 
tions. 

Conclusions:- In summary we have developed a scheme 
for realizing the first single-walker QW in the laboratory, 
with the ion's electronic degree of freedom serving as the 



two-state coin and the motion as the walker's degree of 
freedom. In contrast to current approaches to developing 
QW implementations, which would realize QWs on cir- 
cles in phase space d, Q , our approach yields a RW-QW 
transition in position space. In other words, the walker is 
truly spreading out over unbounded position space rather 
than being folded back on itself. Our approach is true to 
the spirit of RWs over unbounded domains but required 
important innovations taking this idea well beyond TM's 
first concept for the QW in an ion trap 

Although the walk is over position, we show that the 
experimentally accessible phonon number equally reveals 
the RW-QW transition. We have shown that phonon 
number measurement is feasible for dozens of phonons 
by driving the ions at the carrier frequency, then Fourier 
transforming the ground state population to reveal the 
Rabi frequencies ^n,n, hence the phonon number distri- 
bution. This approach is similar to the approach of blue- 
sideband driving but is more effective in revealing 
Pn over wide-ranging phonon number n. In addition, 
we introduce an experimentally controllable phase ran- 
domization procedure that is parameterized by q. The 
RW-QW transition is a key part of any experiment that 
plans to demonstrate QW behavior , yet the ion trap 
dynamics are almost perfectly coherent. Finally we have 
been quite careful in studying LD corrections, which are 
clearly nonnegligible. In conclusion our theory estab- 
lishes a pathway to realizing a many-step QW, and our 
techniques for counting phonons should be useful for gen- 
eral quantum information protocols. 

Note:- Subsequent to submitting this manuscript, 
a trapped-ion three-step coined quantum walk, which 
shows beautifully the difference between the QW and 
RW, has been reported [3. However, their walk is lim- 
ited to three steps to avoid higher LD contributions, and 
they measure motional wavepacket overlap rather than 
position. Our scheme overcomes such limitations. 
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